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ABSTRACT 

Context. We search for persistent extragalactic sources of y-rays with energies above 100 GeV with the Fermi telescope. 

Aims. We construct a systematic survey of the extragalactic y-ray sky at energies above 100 GeV. Such a survey has not been done 

before by the ground-based Cherenkov y-ray telescopes, which have, contrary to Fermi, a narrow field of view. 

Methods. We study a map of arrival directions of the highest energy photons detected by Fermi at Galactic latitudes 1^1 > 10° and 

search for significant point-source-like excesses above the diffuse Galactic and extragalactic y-ray backgrounds. We identify eight 

significant point-source-like excesses in this map. 

Results. Seven of the eight sources are known TeV blazars. The previously unknown source is identified with the galaxy IC 310, 
which is situated in Perseus cluster of galaxies. The source is detected with a significance 6cr above 30 GeV. We identify two possible 
scenarii for y-ray emission from this source. One possibility is that emission originates from the base of relativistic outflow from the 
active nucleus, as in the BL Lacs and FR I type radio galaxies. Otherwise y-ray photons could be produced at the bow shock that is 
formed as a result of the fast motion of the galaxy through the intracluster medium. The two models could be distinguished through 
studying of the y-ray signal variability. 
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Introduction. The number of known astronomical sources 
of very-high-energy (VHE) y-rays grew ten-fold over the last 
five year^. A large part of the newly discovered sources lie 
in the Galaxy and were revealed via a syste matic scan of the 
inner Galactic Plane by the HESS telescope ( lAharonian et al.L 
i2005a, 2006). The HESS survey has covered an area 0.1 sr in 
a strip |/| < 30°, \b\ < 3°. This covers less than 1% of the sky. 
Surveys of larger regions on the VHE y-ray sky with the existing 
ground based Cherenkov y-ray telescopes are difficult because 
the size of the field of view is too narrow (5° for HESS, 3.5° 
VERITAS telescopes and 3° for MAGIC telescope). A previ- 
ous survey of the northern hemisphere using the Cherenkov tele- 
scope Whipple has resulted only in derivat ion of upper limits o n 
the flux of persistent VHE y-ray sour ces ( Weekes et al.L 1979I) . 
The wide field of view MILA GRO dAtkins et all |2004 and 
Tibet dAmenomori et al~ I I2OO5I) arrays have produced a system- 
atic survey of the VHE y-ray sky. However, the energy threshold 
of the air shower arrays like MILAGRO and TIBET is rather 
high (in the multi-TeV band) so that only sources with spectra 
extending well above 1 TeV could be detected. 

Contrary to the ground-based Cherenkov y-ray telescopes, 
Fermi has a wide field-of-view and continuously surveys the 
whole sky on a timescale of 3.2 hr. Over the first year of opera- 
tion Fermi has detected some 1.5 x 10^ Galactic an d extragalac- 
tic so urces of y-rays with energies above 1 GeV dAbdo et al.L 
I2009h . The smaller collection area of Fermi (~ 1 m-, compared 
to ~ 10^ m~ for the ground-based y-ray telescopes) prevents an 



' See e.g. |http : //http : //tevcat . uchicago . edu/ 1 , 
http : //www . mppmu . mpg . de/ rwagner/ sources/ 



extension of the all-sky monitoring with Fermi to the VHE y-ray 
band. 

However, the collection area of Fermi is still sufficient 
for detecting the brightest y-ray sources at the energies above 
100 GeV. The power of the all-sky monitoring capabilities of 
Fermi at the highest energies was clearly demonstrated by dis- 
coveries of new VHE y-ray sources motivated by Fer mi detec- 
tions of these sources above 10 GeV (IOng ll2009Ll20T(ih . The all- 
sky survey capabilities of space-based y-ray telescope EGRET 
at the energies above 10 GeV we re used for the se arch of new 
VHE y -ray blazars bv lDingus & Be rtsch (2001); G orbunov et al.l 
(I2005h via cross-correlation of arrival directions of highest en- 
ergy EGRET photons with positions of known sources. 

Below we use Fermi data to produce a survey of extragalactic 
sky at the energies above 100 GeV, i.e. in the energy range acces- 
sible for the ground-based y-ray telescopes. We find that most of 
the sources visible with Fermi at the energies above 100 GeV are 
known TeV blazars. The only source which has not previously 
been reported as a VHE y-ray so urce turns out to be IC 310 , 
which is a head-tail radio galaxy dSiibring & de Bru^. Il998h 
with possibly a BL Lac type nucleus dRector et al.. 1999). 

Two radio galaxies have previously been reported to 
be t he s ources of y-rays with energies above 100 GeV: 
M87 dAhar onian et al, 2005b; Albert et al., 2008; A cciari et all 
l2009h and Cen A dAharonian et al.Ll2009l) . These two sources are 
the two closest FanarofF-Riley type I (FR I) radio galaxies. The 
FR I rad io galaxies form the "par ent" population of BL Lac type 
blazars dUrrv & Padovani[ll995h . They are expected to be weak 
VHE y-ray emitters, because the y-ray flux from these sources is 
not boosted by the relativistic Doppler effect. In this respect it is 
not surprising that only the two nearest FR I radio galaxies have 
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Fig. 1. Fermi count maps of the sky region around the center of Perseus cluster in 1-10 GeV (left), 10-100 GeV (center) and 100-300 
GeV (right) energy bands. The position of Fermi source IGFL J03 19.7+41 30 is marked by the ellipse corresponding the uncertainty 
of the source position (95% confidence level). The green circle of the radius 0.1° marks the position of additional source detected at 
the energies above 100 GeV. The color scale is linear, from (black) to the maximal number of counts (white). 

been seen in the VHE y-ray band so far. Both Cen A and M87 
are too weak to be detected at 100 GeV in the 1 .5 yr exposure of 
Fermi. 

IC 310 is situated in Perseus galaxy cluster at the distance of 
80 Mpc, which is a factor of 22 and 5 larger than the distances 
of Cen A and M 87, respectively. IC 310 is, therefore, by 1-2 or- 
ders of magnitude more luminous than that of Cen A and M87. 
Besides, IC 310 is not classified as a FR I type radio galaxy. 
Instea d, it is a head-tail radio galaxy (Siibring & de Bruvn, 
Il998l) . the type of galaxies usually found in galaxy clusters. It 
possesses an extended "tail" of the angular size ~ 15' aligned 
along the direction connecting the center of Perseus cluster and 
the galaxy. This tail is believed to be produced as a result of 
the fast motion of the galaxy through the intracluster medium 
dSarazinl Il988h . The detection of IC 310 in the VHE band 
looks less surprising if the hypoth esis that the source hosts a 
low-luminosity BL Lac is adopted jRector et al.L Il999l) . In this 
case the mechanisms of VHE y-ray emission from the source 
could be the same as in the majority of extragalactic VHE y-ray 
sources. 

Survey of extragalactic sky above 100 GeV with Fermi. We 
used the data of Large Area Telescope (LAT) on board the Fermi 
satellite to generate an all-sky map of arrival directions of y-rays 
with energies above 100 GeV. The entire list of photons detected 
by Fermi in this energy range included 4145 events, collected 
over the period from 2008 August 4 till 2010 February 14. A 
substantial fraction of the events is produced by emission from 
diffuse and point sources along the plane of the Galaxy. In our 
analysis of extragalactic sources we exclude the region \b\ < 10° 
from the point source search. This leaves 2603 photons available 
for the analysis. 

To identify significant point sources we first selected clus- 
ters of three and more events within circles of the radius 0.1°, 
corresponding to the 68% containment radius of the LAT point 
spread function at energies above 100 GeV0. For each cluster 
we calculated the probability to find the observed number of 
events within the angular resolution circle by chance. To do this, 
we calculated the overall number A^io of events within a circle 
of 10 degrees radius around the cluster. The chance probabil- 



Table 1. Catalog of y-ray sources detected by LAT at energies 
above 100 GeV at \b\ > 10°. 



ity to find a second event within a 0.1° around the first event is 
p = (l-cos(0.1°))/(l-cos(10°)) ^ lO"! The chance probability 
to find a cluster of at least Kq,\ events within a circle of the radius 



-1 P' 



0.1°isthenP(A?io,/:o.i) = i:A°:^, 
The chance probability Pioo-300 



^(Mo-l)!/(A^io-^-l)!^! 
of finding clusters of Kq,\ > 



^ |http : //www-glast . slac . Stanford. edu/| 
sof tware/IS/glast_lat_perf ormance . htm 



3 events in the 100-300 GeV energy band is less than 10" ev- 
erywhere in the extragalactic sky. The list of /STo.i ^ 3 excesses 
found at Galactic latitudes \b\ > 10° is given in Table [1] In the 
same table we also give the number of y-rays within the circle 
of the radius 0.2°, which corresponds to the 95% containment 
radius of LAT PSF above 100 GeV. In principle, the estimate of 
the chance probability of finding a cluster of > 3 photons in the 
background events could be improved if known Fermi sources 
are removed from the background. However, the only < 10^ 
photons with energies abo ve 100 GeV could b e associated to 
the known Fermi sources jNeronov et al.L 1201 Oh . which means 
that the correction introduced by "masking" the known sources 
in the background estimates would be small. 

All but one detected clusters of Kq i > 3 events listed in 
Table [1] could be identified with the known TeV blazars. The 
only source that does not belong to the blazar class and was not 
previously detected in the VHE energy band is the source #2 at 
RA = 49.14 + 0.08, DEC = 41.30 + 0.06. 

VHE y-ray emission from IC 310. The excess at the position 
of source #2 has three photons at energies above 100 GeV. The 
chance probability that such an excess is due to a chance coinci- 
dence is ruled out with probabiHty (1 - Pioo-300) - 99.9994%, 
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which corresponds to a 4.5cr significance of the source detection 
above 100 GeV. 

Fermi images of the sky region around source #2 are shown 
in Fig. [T] The new source is situated approximately 0.6° to the 
southwest from the bright Fermi source IFGL J0319. 7+4130, 
which is identified with the radio galaxy NGC 1275 in the cen- 
ter of Perseus galaxy cluster The spectrum of the new source is 
harder than the spectrum of NGC 1275, so that emission from 
NGC 1275 dominates in the 1-10 GeV energy band, while emis- 
sion from the new source dominates the signal above 100 GeV. 

At energies above 30 GeV, 95% of events from a point source 
are contained in a circle of the radius 0.3°. This means that 
events from the bright source IFGL J0319. 7+4130 do not affect 
the signal from the newly detected source 0.6° away. Taking this 
into account, we included events with energies above 30 GeV in 
the analysis. We found two more events from the source within 
the 0.1° search circle in the 30-100 GeV band. To calculate the 
probability P30-100 of finding two additional events at the source 
position one can use modification of the expression for the 100- 
300 GeV band with a substitution Kq i - 1 Ko,[ and A^io - 1 ^ 
A^io. Doing this we find P30-100 = 1.1 x 10"'^. The probability 
to find > 2 additional photons in the 30-100 GeV to a cluster 
of > 3 photons in the 100-300 GeV band is Pioo-soo^'ao-ioo - 
6.6 X 10^'" (which corresponds to 6cr for Gaussian statistics). 
Equivalently, the combined probability to find a 3.9 and 4.5cr 
excess in the two bands is ^ 2 x 10"^, which corresponds to the 
significance of source detection at the level of 5.6cr (integrating 
the probability density outside constant likelihood contour for 
Gaussian statistics). 

In order to verify the detection of the source with Fermi, we 
performed a standard Fermi analysi^ for the sky region con- 
taining IC 310. In this analysis we used photons with energies 
between 10 GeV and 300 GeV, collected from a circular area 
with a radius of 5 degrees around the position of IC 310. We 
included all sources mentioned in the Fermi first year catalog 
dAbdo et al.Ll2009l) for this region of the sky. All the sources had 
freely varying normalizations and spectral indices. We assumed 
that the spectra of all the sources included in the likelihood anal- 
ysis had power-law shape. 

The likelihood analysis has resulted in the det ection of IC 
310 with the Test Statistic (TS) (iMattox et al.Lll996l) value of 53, 
which corresponds to approximately 7cr source detection signif- 
icance and is compatible with the 6cr source detection signif- 
icance found from the direct photon counting above 30 GeV. 
Figure |2] shows the map of TS values generated with gtlike tool 
by varying the position of the source on the grid of 12 by 12 
positions with a 0.1 degree step. The uncertainty of the source 
position found from the TS map is also compatible with the one 
found from the direct photon counting. 

Figure |3] shows a comparison of the Fermi image above 
100 GeV with the images of the same region of the sky in 
X-ray (ROSAT a ll sky surveyQ) and radio (WENSS survey 
dRen Relink et al.L Il997h ) bands obtained through the SkyView 
interfac^ The radio and X-ray source at the position of the new 
y-ray source is the galaxy IC 310, which is a part of Perseus 
cluster There are no other bright radio or X-ray sources within 
the 0.1° degree circle around the position of the y-ray source. 
This provides an unambiguous identification of the source with 
IC310. 



' http : //f ermi . gsf c. na sa . gov/ssc/data/ | 

analysis/scitools/ 




Fig. 2. Map of the TS values for an additional source near NGC 
1275 at energies above 10 GeV. The white contour corresponds 
to the 95% uncertainty of the source position. 



In order to estimate the flux from the source in the 30-300 
GeV energy range, we calculated the exposure of Fermi/LAT in 
the direction of the source in this energy band with the gtex- 
posure tool. Collecting the photons from the 95% containment 
circles of the energy-dependent p oint spread function for the 
front and back-converted photons (lAbdo et al l l2009ah . we find 
the estimates of the source flux in two energy bins, 30-100 GeV 



and 100-300 GeV are 



100 



2.TJ.-1 X 10"" erg/(cm^s) and 



' 100-300 



= 1.4 



+1.1 

0.6 



X 10" erg/(cm s), respectively. These esti- 
mates are shown in Fig.|4]together with the multiwavelength data 
on IC 310 collected from the NASA Extragalactic Databas^ 

At energies below 30 GeV, the signal from NGC 1275 could 
give a non-zero number of photons at the position of the new 
source. In order to estimate the background at the source posi- 
tion we adopted the following procedure. We have chosen three 
"background" circular regions of the radius 0.3° situated at the 
same angular distance from NGC 1275 as the new source, but 
at the position angles 90°, 180° and 270°, with respect to the 
direction from NGC 1275 toward the new source. The signal 
of NGC 1275 should give approximately the same number of 
counts in the three "background" regions and in the "source" re- 
gion, which we choose to be a circle of the radius 0.3° around 
the source position. We found that the signal in the source circle 
is compatible at < 3cr level with the signal in the background 
regions in the energy bins 1-3 GeV, 3-10 GeV and 10-30 GeV, 
see Fig.|4] 

Discussion. Known VHE y-ray loud Active Galactic Nuclei 
(AGN) are divided into several classes. Most of the sources are 
BL Lacs, which are relativistically beamed versions of Fanaroff- 
Riley type 1 radio galaxies. Two of the detected sources are 
the nearest FR 1 radio galaxies themselves (M87 and Cen A). 
One source, 3C 279, belongs to the Flat Spectrum Radio Quasar 
(FSRQ) class. Unless IC 310 proves to be a weak BL Lac type 
object, the detection of VHE y-ray emission from a head-tail 
radio galaxy provides a new class of extragalactic VHE y-ray 
sources. If the inte rpretation of the source as a weak BL Lac 
(iRector et al.Lll999l) is adopted, the source still has unusual prop- 
erties due to its unusual for BL Lacs radio mor phology and h ost 
galaxy (lenticular, instead of a giant elliptical) (iNilso 



^ and ho 
nL lT97l . 



[http : //www . x ray.mpe.mpg.de/cgi-bin/rosat/rosat-survey 
' http://skyview.gsfc.nasa.gov/ 



In FR I radio galaxies and BL Lacs the VHE y-ray emission 
is most probably produced in the innermost part of the jet that is 
1- 



|http : //nedwww . ipac . caltech . eduT] 
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Fig. 3. Image of Perseus galaxy cluster in radio band from WENSS sky survey (left), in the X-ray band from the ROSAT all sky 
survey (middle) and the Fermi image in the 100-300 GeV energy band (right). The green contours on the right panel correspond to 
the radio image from the left panel. 
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Fig. 4. Spectral energy distribution of IC 310. Grey points are 
taken from NED. Black points are Fermi measurements. 



ejected by the supermassive black hole. For IC 3 10, it is not clear 
a priori if the y-ray emission is powered by the same mecha- 
nism. An alternative possibility is that y-rays are produced at the 
bow shock formed in interaction of relativistic outflow from the 
fast moving galaxy with the intracluster medium. In this respect, 
the IC 310 - Perseus cluster system might be similar to a much 
smaller scale PSR B 1259-63 system, in which y-ray emission is 
produced at the bow shock formed in interaction of relativistic 
outflow fro m a pulsar moving thro ugh a dense wind of a com- 
panion star dTavani & Aronsl 1997h . 

Angular resolution of y-ray telescopes is only marginally 
sufficient to resolve the bow shock surface in the IC 310 - 
Perseus cluster system. The angular length of the "tail" of the 
bow shock visible in the radio band is ~ 15'. The uncertainty 
of the Fermi source position is ~ 4', which is smaller than the 
length of the tail. The Fermi data indicate that most of the y-ray 
emission is produced in the "head" part of the source (see Fig. [3]). 
At the same time, the Fermi angular resolution is not sufficient to 
distinguish between emission from the "head" of the bow shock 
and the emission from the base of the jet near the supermassive 
black hole, which powers the source activity. 

The crucial test, which would allow to distinguish between 
the two mechanisms would be the (non) detection of variability 
of y-ray emission from the source. Indeed, in the BL Lac type 



models the y-ray emission is expected to be variable at different 
timescales, down to the timescale of light-crossing of the central 
supermassive black hole. On the other hand, if the observed y- 
ray emission is produced at the bow-shaped contract surface be- 
tween the AGN outflow and the intracluster medium, the y-ray 
source has ~kpc scale size. This means that the y-ray emission 
could not be variable on timescales much shorter than ~ 10^ yr 

Variability of the y-ray signal from IC 3 10 could not be stud- 
ied with Fermi, which has detected only five y-rays from the 
source at energies above 30 GeV. We verified that the five de- 
tected events did not come within a narrow time window, which 
would indicate the possibility of a strong flare from the source. 

The presence or absence of variability of the VHE y-ray 
emission from IC 310 could be readily verified in observations 
with ground based y-ray telescopes. A previous observation of 
the region around Perseus cluster with Whi pple telescope has re - 
sulted in an upper limit one the source flux (iPerkins et al.Ll2006l) . 
However, this upper limit is comparable to the Fermi measure- 
ment of the source flux, so that no conclusion about the pres- 
ence or absence of long-term variability of the source could be 
drawn from the comparison of Whipple and Fermi observations. 
It is clear that observations of the source with the new genera- 
tion of ground-based y-ray telescopes VERITAS or MAGIC (the 
source is in the northern hemisphere, not easily accessible for 
HESS) would give much higher signal statistic at energies above 
100 GeV, so that the hypothesis of the flux variability could be 
easily testecfl 

From Fig.[3]is is immediately clear that not every head-tail 
radio galaxy in the Perseus cluster emits in the VHE y-ray band 
at the Fermi sensitivity level. The image of the cluster, shown in 
this figure, includes another prototypical head-tail galaxy, NGC 
1265. This source is clearly identified in the radio band, but, con- 
trary to IC 310, does not show significant X-ray and VHE y-ray 
emission. A comparison of the physical parameters of IC 310 
and NGC 1265 (e.g. velocity through the intracluster medium, 
overall power of relativistic outflow etc.) could help to clarify 
the conditions under which particle acceleration and VHE y-ray 
emission in this type of sources occurs. 



' Following the discovery of the source above 100 GeV with Fermi, 
the detect ion in the same energy band was confirmed by the MAGIC 



telescope l lMariottill20Toh . 
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